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ABSTRACT 



Luminous Blue Variables (LBVs) oftentimes show double-troughed absorption in their strong Ha lines, which are as yet not under- 
stood. Intriguingly, the feature has also been seen in the interacting supernova SN 2005gj, which was for this reason suggested to 
have an LBV progenitor. Our aims are to understand the double-troughed absorption feature in LBVs and investigate whether this 
phenomenon is related to wind variability. To this purpose, we perform time-dependent radiative transfer modeling using cmfgen. 
We find that abrupt changes in the wind-terminal velocity - as expected from the bi-stability jump - are required to explain the 
double-troughed absorption profiles in LBVs. This strengthens scenarios that discuss the link between LBVs and SNe utilizing the 
progenitor's wind variability resulting from the bi-stability jump. We also discuss why the presence of double-troughed P-Cygni com- 
ponents may become an efficient tool to detect extra-galactic LBVs and how to analyze their mass-loss history on the basis of just one 
single epoch of spectral observations. 
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Luminous Blue Variables (LBVs) are unstable massive stars 
whose evolutionary state is under heavy debate. While they have 
been thought to be in a short, transitory phase between an O - 
type and Wolf-Rayet (WR) star (Humphrey s & Davidson! 1994) . 
recent observations suggest that LBVs might actu ally be the final 
stage before the star exp l odes as a tvpe-II SN (e. g. | Kotak & Vinkl 
l2006tlSmifh etai]l2007t iGal- Yam & Leonardll2009h . 

LBVs are characterized by severe mass loss (~ 10 6 - 
10~ 3 M yr~'), slow winds (terminal speeds ~ 50 - 500 kms -1 ), 
photometric and spectroscopic S Doradus-type variations on 
timescales of years, and are often surrou nded by circumstel- 
lar nebulae (Humphreys & Davidson|[r994l) . During the S Dor- 
type variations, the stellar radius changes by up to a factor 
of 10, while the effective temperature changes in the range 
8 000 - 25 000 K dvan Gendered 12001). There is evidence for 
variability in the bolometric luminosity of up to a factor of two 
in the few LBVs that ha ve been analyzed in detail with mo dern 
radiative transfer codes dGroh et al.ll2009h IClark eTal"]|2009l) . 

During the S-Dor variatio ns, LBVs may eventua l ly cross the 
so-ca lled bi-stability limit dVink & de Koterl 120021: iGroh et alJ 
1201 lb . w hich occurs in line-driven winds at around T e ff ~ 
21000^ dLamers et alJ [l99l IVink & de Koterl [2002b. Model 
predictions by Pauldrach & Puis ( 1990) suggested that wind bi- 
stability plays an important role in a shell-ejection scenario of 
PCygni, in the sense that small variations (~ 5%) in L* or 
R+ cause sig nifica nt changes in the wind properties (see also 
iNaiarro et alii 1 997t iHillier et al.lll998l) . 

Today, bi-stability refers to mass loss and wind-ve locity vari- 
ations resulting from changes in the Fe ionization dVink et alJ 
Il999h . Given that the defining characteristic of LBVs are their 
S-Dor variations ((Humphreys & Davidsoni n"994). there is an ex- 
pectation that LBVs could change their mass-loss properties 



due to the bi-stability whilst they trave l back and forth across 
the upper Her tzsprung -Russell diagram d Vink & de Koterl l2002; 
lLamersI 19971) . Indeed. lGroh et alJ (1201 ll) confirmed that the pro- 
totype LBV AG Carinae crosses the bi-stability limit, but were 
unable to identify whether a gradual change or a jump in the 
wind properties occurs. 

It is also the property of having variable mass loss and 
wind terminal speed that first revealed LBVs as candidate 
supernova (SN) progenitors (Kotak & Vink 2006). This as- 
sertion was strengthened by the detection of double-troughed 
PCygni absorption components in the narrow Ha por tion of 
the interacting supernova SN2005gj dTrundle et al. 2008). Such 
multiple-absorption components have been known, for several 
decades, to exist in Ha s pectra of confirmed L BVs (F ig.[p e. g., 
IStahl et al J fl983l l200ll [2001 iLeithereret aflfl994h . The phe- 
nomenon might be the result of s hells ejected by the LBV, or 
caused by radiative transfer effects dHillierll 19921) . or due to vari- 
able winds, as would be expected from the physics of the bi- 
stability jump. 

The main aim of this Letter is to explain the origin of 
the multiple P-Cygni absorption components seen in LBVs. 
Explaining this phenomenon is particularly relevant as it is 
as yet far from proven that LBVs really are SN progenitors. 
A number of hints in this direction have be en suggested, but 
they all have their specific issues (see e.g. iDwarkadasI 1201 ll 
and references therein). Most significantly, stellar evolution the- 
ory does not predict the cores of LBV stars to be advanced 
enough to be anywhere clos e to core-collapse (Lange ret al.l 
ll994tlMevnet & Maede il2003l) . It is thus important to scrutinize 
all suggestions involving a direct link between LBVs and SNe. 
Before being able to attack the line profiles seen in SN2005gj, 
which involve the complexities inherent to both progenitor mass 
loss and the subsequent SNe, it seems logical to first attempt 
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Fig. 1. Observed spectrum o f the LBV AG Ca r around Ha ob- 
tained in 1997 Feb 28, from lStahl etail (120011) . The dashed line 
shows a magnified view (factor of 6) of the double-troughed P- 
Cygni absorption. 



to explain the more widespread existence of multiple absorption 
components in LBVs themselves. 

2. Radiative transfer modeling method 

We employ a modified version of the spherical symmetric, 
fully line-blanketed, n on-LTE radiative transfer code cmfgen 
(Hillier & Millen fl998l) to analyze the effects of time-dependent 
wind parameters on the stellar spectrum. Each model is speci- 
fied by the hydrostatic core radius R+, luminosity L*, mass-loss 
rate M, wind volume filling factor /, wind terminal velocity v^,, 
and chemical abundances Z, of the included species. We com- 
pute time-dependent models with the following stellar and wind 
parameters: T e s in the range of 13000 to 25000 K, M in the range 
of 0.2 to 2 x 10~ 5 M yr-\ M = 60 M , L* = 10 6 L Q , / = 0.25, 
and chemical abundances typical of LBVs, with a He mass frac- 
tion of Y=0.62, N/N =7, C/C G = 0.08, 0/O = 0.025, and solar 
abundance for the iron group elements. These parameters are 
roughly characteristic for AG Car. 

Here we briefly describe our implementation of time- 
d ependent density an d velocity structures. The reader is referred 
to lGroh et al.1 (120081) and Groh & Hillier (201 1, in prep) for fur- 
ther details. At the moment, no time dependence is accounted 
for in the solution of the radiative transfer and rate equations, 
since the recombination timescale is much shorter than the flow 
timescale for the wind regions studied here. As the wind hydro- 
dynamics is currently not yet solved for, in order to account for 
time-dependent outflows, we allow for arbitrary v(r) and p(r) 
stratifications as inputs into cmfgen. We consider cases when the 
star crosses the bi-stability limit, evolving between two epochs: 
the initial epoch 1, when T e s = 23 000 K and a fast wind is 
present with Voo.t, and epoch 2, when T e g decreases and a slow 
wind is present, according to Voo,2 < Voo.i ■ For simplicity, we keep 
M fixed, so that p(r) is computed following the equation of mass 
conservation. 
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Fig. 2. (a): wind velocity law for a time-dependent model at Af = 
120 days with a gradual change in v M (t) (black solid line) and 
with a abrupt change in Voo(f) (red dashed line), in comparison to 
a steady-state model (blue dash dotted line), (b): corresponding 
Ha line profiles. 

We explore two scenarios as to whether the double P-Cygni 
profiles can be caused by the time variability of wind parameters: 
a gradual change in Voo, or an abrupt instantaneous variability. 

For models with gradual changes in Voo(f), for simplicity we 
assume that, starting at t = 0, v^o decreases linearly with time. 
For a given time Af after t — 0, we compute the distance traveled 
by material ejected between t — and Af. The result is a velocity 
law that increases linearly with the distance up to a point where 
the initial wind is reached. From that point onwards, v(r) follows 
the standard /?-type, asymptotically reaching Voo.t (Fig. [2^). 

For models with an abrupt change, we assume that at f = 0, 
Vk, switches instantaneously from v«i,i to v Mi 2. To compute a 
model for a given Af, we evaluate the distance traveled by the 
last particles of the initial wind (r\ ), as well as the distance trav- 
eled by the first particles of the final wind (fz). Following a 
suggestion by S.P. Owocki (2010, priv. comm.), we assume that 
the initial wind characterizes the region from r\ up to the outer 
boundary of the computational domain, whilst a jfi-type law char- 
acterized by Voo.i and f!\ is assumed. From the star up to r2, we 
assume that v(r) is given by the velocity law of the final wind 
according to and @2- We further assume that the region in 
between the two winds, i.e. between n and r\, v(r) and p(r) are 
linear, and join the two winds smoothly. We defer to future hy- 
drodynamical computations for a more accurate treatment of this 
region, but we note that initial tests using somewhat different as- 
sumptions may lead to small changes in the resulting P-Cygni 
line profiles. These changes are not anticipated to affect the con- 
clusions reached here. 

For both gradual and abrupt models, we assume that the stel- 
lar wind is at f = in steady-state with v«, t i = 150 kms" 1 
and y6i = 1 and, as r -> oo, evolves to a final state with 
Vco,2 = 70 kms -1 and = 1- The choice of v x values is based 
on LBV spectral data, whilst the fact or of two chang e in v^, is 
adopted from theoretical predictions (IVink et al.lll999h . 



3. Time-dependent P-Cygni line profiles: gradual vs. 
abrupt changes in the wind parameters 

Figure[2|3 presents the P-Cygni line profiles calculated for both a 
gradual and an abrupt variation of the wind parameters as a func- 
tion of time, here computed at Af = 120 days. A gradual change 
in Voo(f) causes a linear variation of the velocity with distance 
(Fig. [2^). At Af = 120 days, the P-Cygni profile deviates only 
slightly from a steady-state model with the same stellar param- 
eters. To test the robustness of our results, we computed models 
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for different At, as well as stellar and wind parameters all with 
gradual changes, but we still could not obtain double P-Cygni 
absorption components. 

On the other hand, we found that an abrupt change in Vco 
causes significant effects in the wind velocity law (Fig. |2^). 
While the velocity law of the outer wind is similar to that of 
the gradually changing as well as steady-state models, the inner 
part of the wind follows a yS-type law asymptotically approach- 
ing Voo 2. The abrupt change in modifies not only the veloc- 
ity and density structures, but also causes the appearance of a 
double-troughed P-Cygni absorption component (Fig. |2j)). The 
double P-Cygni absorption component arises as a result of the 
presence of a small amount of material with velocities in be- 
tween Vco.i and v M ,2, causing a very low Ha line opacity between 
these two velocities. Models computed with a variable M across 
the bi-stability jump show a similar qualitative behavior as the 
ones shown here, implying the abrupt change in v x is the key 
phenomenon responsible for the double-troughed P-Cygni ab- 
sorption components. We are not able to tell whether there is a 
jump in M as well. 

Figure [3] presents the temporal evolution of the Ha ab- 
sorption profile for selected At, showcasing how the absorption 
evolves from the initial fast wind in steady-state (panel a) to the 
final slow wind in steady-state (panel /, red dashed line). We 
find that the P-Cygni absorption component from the fast wind 
becomes weaker, as the slow wind takes over, disappearing at 
around At = 240 days. The P-Cygni absorption component from 
the slow wind first appears rather weakly at around Af = 90 days, 
and increases in strength as a function of time. Note that even 
when no absorption from the fast wind is seen, its weak emission 
still influences the spectral morphology, even at Af = 960 days. 

So when can the double P-Cygni absorption components be 
seen after the abrupt change in Voo ? The duration of the double P- 
Cygni absorption period depends on the characteristic flow time 
(dependent on R^/v^ and /3) of both winds, as well as on the 
location of the line forming region of interest. For the param- 
eter space explored here, Ha forms at a distance from ~2 to 
10 R+, and the double P-Cygni absorption component is seen 
from At around 90 to 180 days, with its maximum strength 
around 120 days. Higher values of j3 increase the flow timescale 
and allow the double P-Cygni absorption components to be ob- 
served during a more extended period of time. 

4. The dependence of the double P-Cygni 
absorption profiles on r eff and M 

An abrupt change in v m is required to reproduce double P-Cygni 
absorption profiles (Fig. [2]). The next relevant question is: for 
which range of stellar and wind parameters does a double P- 
Cygni absorption profile appear in Ha? 

Here we focus on changes in T e g and M, keeping the other 
stellar and wind parameters constant as well as Af=120 days. 
We found that the morphology of the double P-Cygni absorption 
profiles is strongly sensitive to both T e g and M, and that a zoo of 
absorption line profiles are predicted by our models (Fig. [4]). 

We note that an increase in r eff from 13000 K to 18400 K 
causes less absorption by the initial, fast wind, resulting 
in weaker high-velocity absorption components. Models with 
T e ff > 18400 K have H ionized throughout the wind, which 
dramatically reduces the population of the n=2 level related 
to the Ha absorption. As a consequence, these hotter models 
show extremely weak or no P-Cygni absorption, since the time- 
dependent velocity structure is not probed by Ha if H is ionized 
throughout the wind. 
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Fig. 3. Variations of the synthetic Ha line profiles as a function 
of Af for a fixed r eff = 15.5 kK and M = 9 x 10~ 6 M yr _1 . 
The red dashed line in panel (f) corresponds to the steady-state 
profile obtained as f — > oo. 
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Fig. 4. Variations of the synthetic Ha line profiles as a function 
of r eff (in kK) and M (in units of 1 0~ 5 M Q yr _ 1 ), at Af = 1 20 day s . 



In addition, our model predict weak blue-shifted emis- 
sion between the two absorption components. This blue-shifted 
emission arises in the initial, fast wind, from non-core rays. 
According to our models, the intensity of the blue-shifted emis- 
sion increases as T e s increases. 

Perhaps surprisingly, our models predict that, for a given 
and T e ff , double-troughed P-Cygni absorption appears in Ha 
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Fig. 5. Morphology of the Ha P-Cygni absorption component 
in LBVs as a function of M and r e jf, at At = 120 days. 

only for a relatively narrow range of M. We found that there is an 
optimal value of M (M opt ) when the contrast between the two ab- 
sorption components is maximal. The double-troughed P-Cygni 
absorptions are present when M is about ~ +30% around M opt . 
For Tgff < 18000 K, we found four possible regimes of the Ha 
absorption morphology (Fig. [5]), depending on M: 

• M > 1.3M opt : absorption from the fast wind only, since the 
emission from the fast wind from non-core rays fills out the ab- 
sorption from the slow wind; 

• 0.7M opt < M < 1.3M opt : double P-Cygni absorption, i.e. ab- 
sorption from both winds; 

• 0.3M opt < M < 0.7M opt : absorption from the slow wind only 
since, for lower M, the fast wind is not optically-thick enough in 
Ha to produce detectable absorption. 

• M < 0.3M opt : no P-Cygni absorption, since both winds are not 
optically-thick enough in Ha. 

We note that the ranges given above are only indicative and 
valid for Af=120 days, since the relative location of the four dif- 
ferent regions described above varies with At. In reality, the ac- 
tual observational detection of the double P-Cygni absorption 
should also depend on the spectral resolution and the signal-to- 
noise ratio of the observations, amongst other factors. 

5. Discussion 

We demonstrated via time-dependent, radiative transfer model- 
ing with cmfgen that abrupt wind-velocity changes are needed to 
reproduce the observed double P-Cygni absorption components 
seen in a plethora of LBV spectra. The most promising manner 
causing such abrupt c hanges in v„ is by mea ns of the bi-stability 
jump, as discussed in I Vink & de Koten (120021) . 

We reinforce that the modeling performed here is not ex- 
haustive, given the vast parameter space to be explored and 
the model limitations. For instance, the full hydrodynamics and 
time-dependent solution of the radiative transfer and statistical 
equilibrium equations are not yet accounted for. Furthermore, we 
consider an idealized scenario of an instantaneous abrupt change 
in Voo, whereas in reality, the changes could occur over extended 



time periods of days or weeks. Constraints on the timescale 
could be obtained via time-dependent modeling of observational 
data of LBVs that actually cross the bi-stability limit, e.g. AG 
Car. These topics will be the subject of future studies. 

We note that the empirical identification of a luminous mas- 
sive star as a bona-fide member of the LBV class is a challeng- 
ing exercise. This is mainly due to the long timescales (years to 
decades) necessary to monitor the irregular spectral variability. 
As the detection of double P-Cygni absorption components is 
a strong indication that the candidate object is an LBV under- 
going S-Dor type variability cycles - with its abrupt changes in 
Vk, on timescales of years - we suggest that the sheer presence 
of such double P-Cygni absorption components may become an 
efficient, novel way to detect Galactic and extra-galactic LBVs. 

We found that the morphology of the double P-Cygni absorp- 
tion components depends strongly on M and T e ff . This opens up 
the possibility to constrain these quantities without the lengthy 
process of non-LTE radiative transfer modeling: once L* is 
known, for given At, Vco.i, and Vco,2 could be derived empirically. 

If observations at multiple epochs are available, we should 
also be able to trace the LBV mass-loss history using a similar 
modeling strategy as employed here, as the Ha line forms over 
an extended stellar-wind region in LBVs. This will allow us to 
use the double P-Cygni absorption components as a tool to look 
back in time, and to be able to estimate LBV stellar parameters 
from one single epoch. Although in this Letter we focused on 
Ha, motivated by the observational data available, we note that 
our time-dependent models predict similar double-troughed P- 
Cygni absorption in many other spectral lines formed at different 
wind regions, such as most of the Balmer lines, Fe II, N II, He I 
/U0830, among others. Analyses of these lines should lead to a 
much improved understanding of the mass-loss history of LBVs. 

Interestingly, the morphology of the absorption-line pro- 
files seen in Fig. |4] closely resemble t hose found in the inter- 
acting SN 2005gj dTrundle et all 120081) . This strengthens the 
link between LBVs and SNe that based their assertions on 
the b i- stability mechanism during S Dor cycles (Kotak & Vinkl 
l2006t rTrundle et al l 120081 and the LBV/SN link in general 
(ISmifh et all 120071: iGal-Yam & Leonard! l2Q09h . We expect that 
our results will encourage further observational and modeling ef- 
forts to establish under which conditions LBVs explode as SNe. 
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